Spin-selective localization due to intrinsic spin-orbit coupling 
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We study spin-dependent diffusive transport in the presence of a tunable spin-orbit (SO) interac- 
tion in a two-dimensional electron system. The spin precession of an electron in the SO coupling field 
is expressed in terms of a covariant curvature, affecting the quantum interference between different 
electronic trajectories. Controlling this curvature field by modulating the SO coupling strength 
and its gradients by, e.g., electric or elastic means, opens intriguing possibilities for exploring spin- 
selective localization physics. In particular, applying a weak magnetic field allows the control of 
the electron localization independently for two spin directions, with the spin-quantization axis that 
could be "engineered" by appropriate SO interaction gradients. 

PACS numbers: 72.25.Dc,73.20.Fz,71.70.Ej,85.75.-d 



I. INTRODUCTION 

Recent years have seen a growing interest in spin-orbit 
(SO) interactions in metals, semiconductors, and topo- 
logical insulators, with significant advances in the physics 
of spin Hall effects and other phenomena, whereSO cou- 
pling profoundly modifies electronic transport! 1 * 2 * 3 * 4 ! The 
interest in these problems is also fueled by the desire to 
develop spintronic circuits that do not rely on magnetic 
elements. 

Spin-transport properties are in general associated 
with a gauge freedom of position-dependent spin ro- 
tation. This SU(2) gauge symmetry and the associ- 
ated curvature field were recently discussed in the con- 
texts of the spin Hall physics, the Aharonov-Casher ef- 
fect, and other spin-transport phenomenaP Additionally, 
there have been parallel studies on laser-induced non- 
Abelian gauge fields in cold-atom optical lattices P which 
discussed the Hofstadter "moth" spectrum, the Ander- 
son localization transition, and the possibility of realizing 
non-Abelian interferometry. Understanding the struc- 
ture of the gauge-covariant curvature field underlying the 
SO interaction may thus carry fundamental and broad 
implications. 

In this article, we explore the physical consequences 
that arise from topological structure of this SU(2) gauge 
field. As a practical example, we will consider quantum 
transport corrections in a disordered two-dimensional 
(2D) conductor with intrinsic SO coupling. In particular, 
by focusing on the Rashba interaction we demonstrate 
the possibility of a tunable spin-selective localization of 
the charge carriers, by spatially modulating the strength 
of the SO interaction, which can be achieved via the care- 
ful control of electric gates or by elastic strain. 

Our formalism is based on the geometric nature of the 
precession of electrons in the presence of SO fields, which 
will facilitate the discussion of the weak-localization ef- 
fects: The underlying quantum interference between 
closed time-reversed trajectories can be analytically de- 
scribed by a path integral over a non-Abelian gauge field 
or Wilson-loop integral. We find that the effective size 
of these pairs of trajectories provides a natural limit for 



separating Abelian and non-Abelian contributions. A 
compact, Abelian treatment using the Stokes' theorem 
can be applied to smaller loops with a moderate net spin 
precession, while a large spin precession is responsible 
for the destruction of the phase coherence of longer time- 
reversed loops, where the non-Abelian contributions play 
a dominant role. Our discussion is relevant to recent ex- 
perimental studies of weak (anti)localization in 2D het- 
erostructures with strong intrinsic SO coupling, 8 and is 
complementary to the established theoretical framework 
for disordered systemsP by explicitly revealing the topo- 
logical structure and exploring the role of SO gradients. 
In the following, we will start by formally casting the 
SO interaction as a Yang-Mills gauge field and the in- 
terference effects in terms of the gauge-covariant curva- 
ture field, using the generalized Stokes' theorem. This 
SU(2) curvature can be diagonalized in spin space, with 
its eigenvalues defining a spin-dependent fictitious mag- 
netic field. The corresponding spin-quantization axis 
has in-plane and out-of-plane components, that depend 
on both the strength and gradients of the SO coupling. 
Therefore, the SO interaction combined with an external 
magnetic field generates a total effective field that can 
be tuned independently for spins up and down along a 
desirable spin-quantization axis. These ideas can be ap- 
plied towards a semiclassical description of SO effects on 
the weak-localization physics in 2D, providing a simpler 
and physically more intuitive picture than offered by the 
conventional diagrammatic Cooperon perspective. Our 
approach is especially advantageous when dealing with 
system having inhomogeneous SO interactions. 



II. MODEL 

A. Effective SU(2) gauge field 

Consider the Hamiltonian for 2D electrons in the pres- 
ence of both electric and magnetic fields. An SO interac- 
tion arises from the coupling of an electron's spin with its 
orbital motion, which stems from relativistic corrections 
as electron accelerates in the presence of crystalline or 



applied electric field. Although there can be many differ- 
ent forms of the SO interaction, the primary difference 
originates from the type of electric field source. For ex- 
ample, when the source electric field is from an in-plane 
impurity potential, the type of SO interaction is called 
"extrinsic." In contrast, an applied electric field, elastic 
strain, crystalline inversion asymmetry, or 2D confine- 
ment asymmetry perpendicular to the two-dimensional 
electron gas (2DEG) results in an "intrinsic" SO interac- 
tion incapsulated in the following Hamiltonian: 



F=- 1 -(-/7/V 

2 m 



A(r) -eV(r) 



(1) 



where A(r) is an atomically smooth function, m is the 
effective mass, — e is the electron charge, and V(r) is a 
scalar in-plane, quenched impurity potential. Hamilto- 
nian (111) follows from the effective-mass approximation, 
and we retained the general linear in momentum spin- 
orbit interaction, in terms of the 2x2 "vector potential" 
A, or "connection" with noncommuting matrix compo- 
nents. The associated "curvature" is described by the 
covariant tensor fieldP^ 



J" av — du,A v On A, 



2ni 



|A M , A v \ 



(2) 



where /i and v label the two spatial components, and 
we have, for simplicity, neglected Zecman interactions 
(which would produce also the time component of the 
connection field). </> = hc/e is the magnetic flux quan- 
tum. It is easy to verify that a position-dependent SU(2) 
wave-function transformation corresponding to Hamil- 
tonian (Ilk, ip = t/(r)i//, transforms the connection as 
A' = U^A^U — (i(po/2Tr)U'd l j,U, and the associated cur- 
vature simply as T' v — WT^U. There is also the 
usual U(l) symmetry associated with the ordinary (spin- 
diagonal) contribution to the vector potential, which de- 
scribes the magnetic field. 



B. Wilson-loop integral 

The 2x2 single-electron propagator from r^ to r/ along 
a certain spatial contour C in time t, corresponding to the 
Hamiltonian (nj), is given by 



K(r f ,r i ;t)=T c e 



-(2iri/0 o )/ c dr-A(r) 



K(T f ,T ti t). (3) 



Here, K(rf,Ti',t) is the spin-diagonal propagator if we 
set A = and Tq is the path-ordering operator along 
the contour C. If we are interested in the interference 
along different trajectories, the basic building block for 
describing SO coupling effects is provided by the Wilson- 
loop integral: 



W c {v) = T c e 



-(2iH/0 o )£,dr'-A(r') 



(4) 



which, for general spatially inhomogeneous systems, is a 
function of the position r and the closed contour C, start- 
ing and ending at r. Notice that Wq is a purely geometric 



object dependent on the contour but independent of how 
fast the particle moves around it. Furthermore, to clar- 
ify our nomenclature, this object is called a Wilson-loop 
integral rather than simply a Wilson loop, because the 
trace operation has not been performed. 

We shall demonstrate in the next section that the the- 
ory simplifies either in the presence of strong SO spatial 
inhomogeneity, making it effectively Abelian to the lead- 
ing order in the SO strength,™ or in the case of a simple 
geometry, such as a single-loop Aharonov-Casher effect. 4 
The limit of mesoscopic systems small on the scale of the 
spin-precession length is also simple from our point of 
view, since only the leading non-Abclian corrections to 
the geometric spin transformation need to be retained. 
Fortunately, furthermore, we will argue that the physics 
of weak localization also does not require the full non- 
Abelian treatment in many cases of practical interest. 
Besides, the weak localization provides an illuminating 
example, where we can easily relate the Wilson-loop per- 
spective at SO effects with readily measurable quantities. 



III. THE NON-ABELIAN CORRECTIONS 

The Wilson-loop integral Q provides a gauge- 
covariant description of spin precession and interference 
in the presence of an SO interaction. In the case of an 
Abelian theory, such as U(l) clcctromagnetism, we can 
employ the regular Stokes' theorem for evaluating Eq. (El 
as follows: 



W e (r) 



-27T10/0O 



(Abelian) 



(5) 



where <j> = J dS ■ (V x A) is the magnetic flux integrated 
over the oriented surface area dS. Note that the Stoke 's 
theorem, Eq. (I5J) , provides the manifestly gauge-invariant 
form for We, expressed in terms of the magnetic field 
B = V x A. On the other hand, the non-Abelian Stokes' 
theorem is required for the SU(2) theory, which relates 
the curvature given by Eq. (pi to the exponentiated loop 
integral given by Eq.Q, resulting in the following surface 
integral: 



Wc(r) = T s e -O™/0o)/ s dr'^Adr'^^ir') 



(6) 



Here, T s is the surface-ordering operator^ and T^ is the 
path-dependent curvature, defined by 



T ilv (v')=k(v,v')f ilv (v')K(v\v) 



(7) 



Clearly the mathematical structure for non-Abelian fields 
requires a more elaborate approach than the familiar 
Abelian version, and one has to contend with an ob- 
ject that is more complex than a simple magnetic flux. 
A more general 2D theory would require evaluating 
Eq. (f6J) over many possible closed trajectories, which is a 
formidable task without a simple representation for the 
surface integral on the right-hand side of Eq. (fsl) . There- 
fore, as a first step, let us examine some asymptotic limits 



within this framework. For simplicity, we will take spe- 
cific case of the Rashba form, corresponding to 



the fictitious magnetic field, 3 



A(r) = A(r)z x <r . 



(8) 



where r lies in the xy plane, z is a normal unit vector and 
cr is a vector of the Pauli matrices, which are generators 
of the SU(2) group. A schematic of our model is shown 
in Fig. [II 



A. Homogeneous SO coupling 

As the SO coupling strength A(r) becomes more ho- 
mogeneous, the derivative components of Eq. (pi) van- 
ish and the curvature field is determined entirely by the 
commutator: T\ 2 oc i[Ai,A2] oc & z . Consider the in- 
tegration in the exponent of Eq. (H|. For simplicity, 
we can take the closed contour C to be a unit square 
loop in the xy plane. For a uniform square, connecting 
the points (0,0) -> (1,0) -> (1,1) -» (0,1) -» (0,0), 
Wc — e lAy e lA ' c e~ lAv e~ lA * (absorbing the 27r/0 o fac- 
tor by the redefinition of the connection A). Since the 
Cartesian components of A are noncommutative, the 
concatenated exponents must be expanded in the Baker- 
Campbell-Hausdorff series, which to the leading nontriv- 
ial order in SO strength becomes: 



We = e iAy e iA * e~ iAy e iA * 

_ iA x +iA v +[A lc ,Ay]/2+- e -iA m -iA v +[A x ,Ay]/2+- 



(9) 



In this approximation, Wc(r) can be rewritten in the 
form of the Abelian Stokes' theorem as (restoring the 
2ir/cf>Q factor): 



Wr 



D — 2Tricj>/(f>o 



(10) 



in terms of the "plain" flux of the field T through the loop 
in 2D, without any surface ordering or path dependence, 



dr'» A dr lv T^{v') = / dS T 12 {v') . (11) 



1 



For the Rashba model (pi), this becomes^ 

4> = -{An\ 2 S/fa)cr z = -7r<Mr/Z so ) 2 ,7 z , (12) 

where l so — 4>q/2X is the spin-precession length. The 
relative correction to the plain flux in the exponent of 



Eq. ( 10 1 scales as r/l so (irrespective of the detailed shape 
of the loop) , so that the approximation requires that r <C 

'so • 



B. Strongly Inhomogeneous SO coupling 

For strong spatial variations in A(r), the first term of 
Eq. (pi dominates, and such a variation of A(r) results in 



B = V x A = (& ■ VA) z, 



(13) 



directed along the z axis. With respect to the spin space, 
the fictitious field B has opposite signs for spins up and 
down along the gradient VA. In Ref. 3, the problem 
of the "boundary spin Hall effect" was discussed where 
A(r) changes abruptly from some constant value Ao to 
zero along the edge of a 2D Hall bar. This problem is of 
interest as a candidate model for a lateral Hall contact for 
spin injection by a Rashba system.QU If the SO strength 
A variation is abrupt on the scale of l so , the fictitious field 
( |13[ ) dominates the boundary physics, resulting in a spin- 
dependent (ordinary) Hall effect. The noncommutative 
contribution to the gauge field (pi) in this case is small. 

The lengthscale characterizing the SO inhomogeneity 
is given by 



/inh ~ A/|VA| 



(14) 



which has to be shorter than Z so for the gauge struc- 
ture to reduce to the essentially Abelian theory. In such 
strongly inhomogeneous limit, local (on the scale of l so ) 
SO physics is dominated by the fictitious field ( 13 ), while 
the non-Abelian commutator contribution to the curva- 
ture (p| becomes relatively unimportant. Invoking the 
Abelian Stokes' theorem, we thus reduce the problem to 
the approximation (10), where the flux is now given by 
the familiar relation: 



[ dS-B= J dS{a-V\) 



(15) 



By comparing the limits of the homogeneous vs 
strongly inhomogeneous SO interactions of the Rashba 
model, one discovers one conspicuous difference: In the 
former case, the spin quantization axis set by the cur- 
vature "flux" (12) is along the normal z, while in the 



latter it is along the i n-p lane direction determined by 
the gradient VA, Eq. (15 1. 4 In the very special case of 
a combination of the linear Dresselhaus and Rashba SO 
interactions of equal strength, we have an Abelian SU(2) 
field with zero curvature in homogeneous systems. 5 . In 
this trivial limit, the SO coupling can be eliminated by 
an appropriate gauge transformation. 



C. Small-loop limit 

Both the homogeneous and strongly inhomogeneous 
limits con side red above reduced to the same approxima- 
tion, Eq. ( 10 ) , in terms of the 2x2 flux <fi of the curvature 
tensor component J-\ 2 . The difference between the two 



limits was only in which was the dominant contribution 
to the curvature Q2J> : the derivative, <9iA 2 — d 2 Ai, or the 
commutator, [Ai,A2], piece. In particular, in both of 
these special cases, we were not concerned with the sur- 
face ordering required by the non-Abelian Stokes' theo- 
rem, Eq. ([6]), as long as spins precess and interfere on 




^so 



FIG. 1: A schematic of the considered two-dimensional model: A 2DEG residing in the xy plane experiences a position- 
dependent SO coupling parametrized by \(x,y). The dark patch shows a region where the SO strength is enhanced, e.g., by an 
electric gate or mechanical strain. The SU(2) phase factor accumulated by spin precession in the SO field is r epresented by the 
path- and origin-dependent Wilson- loop integral Wc(r), Eq. (Ml. In the small-loop approximation, Eq. (10 1, the Wilson-loop 
integral is directly related to the Yang-Mills curvature "flux" <j> through the loop. The curvature (|2j is dominated by the 
derivative terms, <9 m j4„ — dvA^, in the regions where the SO strength X(x, y) is rapidly varying (e.g., close to the edges of the 
dark patch), while its non-Abelian nature is reflected in the commutator contribution, L4 MlJ 4„], in the regions where the SO 
strength is smooth or constant. The small- loop approximation requires that r <C l ao , where r is the characteristic size of the loop. 
The weak-localization transport corrections near some point r are governed by the interference between counter-propagating 
closed trajectories starting and ending at r. 



lcngthscales less than the characteristic l so . In the op- 
posite limit, the fully non-Abelian topological structure 
should become manifest. 

Let us notice that the integral Q enters naturally in 
the semiclassical path-integral construction of the quan- 
tum transport corrections P^l In the next section, we will 



for diffusion, J = — -DVp, where the ordinary diffusion 
constant in two dimensions is: 



demonstrate that Eq. (10) gives an adequate approxi- 
mation for many cases of interest, in this context. The 
physical reasoning behind this is as follows: The onset 
of the non-Abelian effects as the loop size approaches 
l so also signals the onset of destructive spin interference 
between time-reversed trajectories and thus suppression 
of quantum corrections stemming from the larger loops. 
The singlet Cooperon channel, which is invariant under 
SO precession, provides an exception to this argument. 



IV. TRANSPORT CORRECTIONS IN A 2DEG 

A. Preliminaries 

We start by recalling the path-integral formalism be- 
hind the semiclassical treatment of weak localization and 
the DC transport in disordered systems. Let us consider 
a system of noninteracting electrons moving in a random 
environment. Clearly, any deviations from an ideal crys- 
tal will result in scattering events that will contribute 
to the resistivity. The classical paths of the scattered 
electrons can be represented by Brownian motion. In a 
2DEG, the electrons will move around the material with 
group velocities close to the Fermi velocity vf- Further- 
more, in the limit of a large system size, the particles 
that experience a random walk will also obey Fick's law 



D = v%t/2 = 1 2 /2t , 



(16) 



for a given mean free path I and scattering time r. The 
conductivity a and D are related by the Einstein relation 
a = e 2 N(0)D, where N(0) is the density of states at the 
Fermi level. It should be emphasized that in the semi- 
classical approach, D can be computed explicitly from 
the velocity-velocity correlation function, which can be 
understood as characterizing the time t it takes for an 
electron to forget its initial velocity direction. Since we 
are interested in the contributions to the velocity- velocity 
correlation function that are responsible for quantum cor- 
rections to the DC transport, we focus on the set of the 
Brownian paths the electrons will take and the associated 
probabilities for propagation and return. 

Weak localization is a quantum-interference correction 
to a that results from phase-coherent backscattering. 
Consequently, these corrections follow from a modified 
diffusion model that includes interference effects. Any 
additional phenomena modifying the weak localization 
itself must be corrections affecting time-reversal proper- 
ties of electronic trajectories, which can destroy the phase 
coherence of back-scattered waves. In particular, since 
the singlet Cooperon channel is responsible for the weak 
antilocalization, suppressing the triplet channels may de- 
localize electrons. The most commonly used formalism in 
weak localization calculations is usually taken in the lan- 
guage of impurity averaged propagators and maximally 
crossed Langer-Neal diagrams. However, a semiclassical 
approach is more intuitive and physically appealing, as 
the effect can be understood in a more transparent de- 



scription as a modified diffusion process yielding the fol- 
lowing quantum correction to the classical conductivity 



4e 2 



dtR{t) 



,-t/r 



(17) 



Here, dtR{t) is the return probability associated with 
the interference of all classical paths with their time- 
reversed counterparts within the time interval (£, t + dt). 
The elastic scattering time r serves as the lower limit 
of integration, while the upper limit is effectively set 
by the phase-coherence time t v 3> r, according to the 
factor e.~ t l T f . Furthermore, the explicit averaging over 
the random potential can be carried out by incorporat- 
ing random fluctuations into the classical paths. The 
classical probability of return can be formally written 
as R(t) = f d[r(t')]P t [r(t')], for closed paths in a neigh- 
borhood of the trajectory r(i'). The probability for the 
realization of a (coarse-grained) Boltzmannian path r(t') 
is given by the Wiener measure: 



PtW)] 



,-f*dt'r{t') 2 /AD 



(18) 



In the presence of a SO interaction (111, the interference 
between counter-propagating closed trajectories acquires 
an SU(2) phase correction Q. For a sufficiently weak SO 
interaction, the trajectories are assumed to be otherwise 
unaffected. (We will specify the exact condition below.) 
The interference between the counter-propagating closed 
trajectories is affected by the spin precession as follows?^ 

R(t) = j ' d[v{t')]P t [v{t')]WcW)]wl[v{t')] , (19) 

where f(i') is the time-reversed trajectory. The re- 
turn probability R entering Eq. ( 17 I has to be appro- 



priately spin-averaged: R = Tr[i?]/2. Let us notice that 



W,J[r(£')] = Wc[r(£')] and use the approximation (10) to 
find for the spin return interference: 



WcW)]wl[r(t')] » e-*** 1 * 



(20) 



We will return later to discussing the range of validity of 
this approximation in the context of weak localization. 
Tracing this over spin to get the quantum conductivity 
correction ( 17 1 , we see that the calculation now reduces 



to finding the eigenstates ±(/> of the traceless Hermitian 
matrix <f>, in the case of an SU(2) field. For the more 
general SU(2)xU(l) field (e.g., Rashba SO plus the ordi- 
nary electromagnetic vector potential) , cf> acquires a finite 
trace, with the two eigenstates <j>± becoming 



= </>m ± ' 



>er • n. 



(21) 



Here, </> m is the ordinary magnetic flux through the loop, 
</> is the SO contribution, and n is a unit vector defining 
spin-quantization axis. In general, by not relying upon 
the approximation (10), the same composition is upheld, 



although the Hermitian matrix <fi that determines spin 
prec essio n for a given loop is not given simply by the 
flux (111. 4 Note that we have neglected Zeeman splitting, 



focusing on the case when the magnetic field is oriented 
perpendicular to the plane of the electron motion!^ 

It is well known that weak localization is affected by 
the presence of magnetic fields and SO interactions. Both 
effects modify electronic diffusion and subsequently, the 
phase coherence between time reversed paths. In the 
case of a magnetic field, phase coherence is destroyed 
in a process called "anomalous magnctoconductance." 
Within our Abelian approximation (10), the intrinsic SO 



quantum corrections to the conductivity can be expressed 
in the same language as magnetoconductance due to a 
spin-dependent fictitious magnetic field. However, be- 
cause of the different effect on the singlet and triplet 
Cooperon channels, a sufficiently strong SO interaction 
can reverse the weak localization leading to antilocaliza- 
tion, and the conductivity that would decrease with an 
increasing magnetic field. In this case, the relative spin 
orientation of the interfering electron trajectories is im- 
portant. In particular, in certain special cases of interest, 
it can be possible to localize spin species asymmetrically, 
with the corresponding spin-quantization axis tunable by 
a combination of the strength and the gradients of the 
SO coupling parameter A(r). We call this limit "spin- 
selective localization." 



B. Spin-selective localization 



One intriguing consequence of Eq. (21 1 is the possibil 



ity to delocalizc electrons spin-selectively by an appropri- 
ate combination of an applied magnetic field and an ad- 
justable SO interaction. The relevant spin-quantization 
axis n can, for example, be chosen to be along the in- 
plane gradient of a strongly inhomogeneous SO strength, 
Eq. (15). This provides a pragmatic scenario, since we 



could have large fictitious fields with moderate magni- 
tudes of the SO interaction, while also allowing for the 
in-plane freedom in choosing the spin-quantization axis. 
In practice, however, this hinges on the ability to control 
large gradients of the SO strength with elastic strain or 
electric gates, for instance. 

The quantum co ndu ctivity correction in the presence 
of an effective flux (21 1 is given in 2D for each spin by^ 



.. - ../ L i ^ T * i ..^ 



a =-ip 



2 8tt|0±| t 



1 



2 87T&fc 



(22) 



in terms of the digamma function ijj, in units of e 2 /2ith. 
When \4>±\ <§C <f>o, this approximates to 



-In 



t> 2 (2tt<J)± 
V 3 



(|</»±|«0 O ). (23) 



(j>± entering these equations is given by the flux of the 
effective field B± — B m ± B through the area Z 2 = Dt v . 



B m is the out-of-plane magnetic field, and the fictitious 
field ±B is given by the eigenvalues of the curvature 



^12 



VA - (^o/Z s 2 > 2 



(24) 



assuming for simplicity a constant gradient VA, on the 
scale of l v . Notice that because of the covariance of the 
Yang-Mills curvature J- ', B± must be gauge invariant. In 
the opposite limit of \<p±\ = \B±\l%, ~^> 0o (when, e.g., 



oo), Eq. (22 1 approximates to 



-In 



in\B ± \DT 



(l«£±l»<M 



(25) 



The above limits are obtained easily by using the asymp- 
totic behavior of the digamma function: ip(l/2 + z) = 
\nz + l/2Az 2 + 0{l/z A ) at z -> oo and ^(1/2 + z) = 0{z) 
at z — ► 0. 

It is necessary to rem a rk t hat our scmiclassical treat- 
ment based on Eqs. ( 18 )-( 20 1 requires that \B±\ <C 0o/' 2 
(where I = vft is the mean free path), so that we always 
have (4>o/\4>±\)( T ip/ T ) ^" 1- I n particular, for the field 
B± determined by the SO curvature ( 24 ) , this translates 
into the requirement I <C ( S oi when A is uniform. Note 



that for a fixed and finite SO interaction, Eq. (23 1 cor- 



responds to a dip of the total conductivity as a func- 
tion of the physical magnetic field B m at zero field, since 
(</>% + ^_)/2 = 2 + < 



0, since |<^>+0_| 



B„ 

discussion below leading to Eq. (126 ) 



while Eq. (25) has a peak at 



[See, however, the 



It is now appropriate to discuss the legitimacy of the 
approximation (10 1, in the context of the weak local- 
ization corrections. As mentioned earlier, this approx- 
imation requires that r <C l so . Otherwise, the omitted 
non-Abelian corrections to the Stokes' theorem due to 
spin precession become appreciable. Let us turn off the 
magnetic field (cj> m = 0) and return to the two extreme 
cases discussed in Sees. |III A| and |HI B| Homogeneous A, 
such that the curvature field JF 12 is determined by the 
second term in Eq. (24 1, and large constant gradient 



VA, such that the curvature is dominated by the first 
term. In the former case, the flux through the area l 2 Q is 
4> ~ \ 2 l 2 /<f>Q ~ </>o, which means the non-Abelian correc- 
tions to the preceding weak-localization analysis based 
on Eq. ( 10 1 can become appreciable if l so <C l v . In 



fact, we in general cannot reduce the problem to uncor- 
related propagation of two spin projections, since only 
the triplet subspace of the Cooperon precesses on the 
scale of l so , while the singlet channel is unaffected by 



SO coupling .1121141 g q ( 25 ) corresponding to l so -C I 



pi 



the case of a homogeneous Rashba SO with (f> m = 0, thus 
should only be valid for the triplet channels, for which Z so 
sets the cutoff lengthscale for trajectories contributing to 
the coherent backscattering. Separating the antilocaliza- 
tion singlet contribution, we thus immediately find 



ctq/2 » \n(l v /l) - 3ln(l so /2irl) (l v > l so ) 



(26) 



in agreement with Ref. 14. The factor of 3 accounts 
for the triplet degeneracy. In particular, note the 



crossover from localization to antilocalization, as l so is 
made shorter. In the other extreme of a large spin-orbit 
gradient, we should generally have a similar concern of 
non-Abelian corrections beyond Eq. (10): Although the 



larger gradient increases the curvature field and shrinks 
the relevant lengthscale corresponding to (f> ~ <po, it is 
easy to show the spin precession will remain significant. 
We will not pursue this problem in detail here, but it 
is worthwhile to remark that it can be overcome, e.g., 
by tuning the Rashba parameter A in combination with 
the linear Dresselhaus SO interaction, to reduce the non- 
Abelian contribution to the curvature fieldP' 



C. Intrinsic vs extrinsic SO coupling 

It is also instructive to compare the limit of weak in- 
trinsic SO coupling, Eq. (23 1, to the delocalizing correc- 



tion in the presence of a weak extrinsic (random) SO 
interaction due to quenched disorder!^] 



ctq/2 « - ln(r v /r) + 2t v /t S( 



(27) 



where r" 1 is the extrinsic SO scattering rate due to impu- 
rities. In the case of a homogeneous Rashba SO coupling 
(|SJ> , the D'yakonov-Perel (DP) spin-relaxation rate 1 is 
Tqp tx A 2 r, in the relevant regime of ( < ! so . The lead- 



ing SO correction to the localization ( 23 ) is thus propor- 
tional to Tq P , i.e., the square of the DP spin-relaxation 
rate, since the curvature defining the effective flux <p± is 
proportional to A 2 , for a uniform A. 

We can easily understand the different localization de- 
pendence on the spin-relaxation rates, r" 1 and r^p, in 
the two cases, by examining how the SU(2) phase prefac- 
tor in the propagator d3J) depends on the strength of the 
SO coupling for closed loops. The uncorrelated quenched 
SO disorder leads to a memoryless (Markovian) spin pre- 
cession, which is similar for both open and closed trajec- 
tories. However, the intrinsic SO coupling (Tsl> combined 
with ordinary momentum scattering results in a DP spin- 
precessional random walk, which is qualitatively different 
for open and closed paths. In the case of open trajecto- 
ries responsible for DP spin relaxation, this random walk 
is Markovian. However, the net spin precession is re- 
duced by closing a trajectory, because (A-dr/dt)\ oop = 0, 
hence the higher-order scaling of the delocalizing correc- 
tion with the spin-relaxation rate. Note that this vanish- 
ing of the average field driving spin precession along the 
closed loops is exact only for a homogeneous A. Making 
A inhomogeneous would enhance the delocalizing correc- 
tions, since now (A • dr/dt)\ oop ^ 0. For strongly inho- 
mogeneous A, the curvature (24 1 scales linearly with the 



SO strength A, and the delocalizing correction linearly 
with the DP spin-relaxation rate t^ p . 

This is analogous to the effect of SO interaction on 
weak localization and quantum conductance fluctuations 
in a chaotic quantum dot™ The relevant SO scatter- 
ing rate is substantially reduced due to the geometri- 
cal constraint on particle trajectories through a chaotic 



quantum dot, in the case of a uniform SO vector po- 
tential pj. By introducing SO nonuniformity, however, 
the rate of SO scattering can be augmented to a level 
comparable to that in the bulkPS SO inhomogeneities 
are in fact inevitable even in high-mobility quantum 
wells due to fluctuations in the concentration of remote 
dopant ions. 16 Their important role has been established 
experimentally 17 with regard to spin relaxation in sym- 
metric quantum wells, where the ordinary D'yakonov- 
Perel mechanism is inactive. The nonexponential spin 
relaxation due to random SO field in the presence of mag- 
netic field has also been recently proposed,^ noting, in 
particular, qualitative difference between the open and 
closed trajectory contributions. These interesting effects 
remain to be explored in the context of quantum correc- 
tions to transport. 



tx A 2 with the spin-quantization axis out of plane, while 
the gradient VA controls the in-plane curvature compo- 
nent. The combination of the two determines the net 
spin-quantization axis n, while the eigenvalues ±B of 
the curvature matrix T\i, which are opposite for the two 
spins, combined with the ordinary magnetic field B m , 
which is the same for the two spins, can in principle re- 
sult in any desirable spin-dependent field B±. One in- 
teresting possibility would be to tune the field B± deep 
in the localized regime such that it vanishes for one spin 
species only, say B^ =0. In this case, only the spin f 
would be delocalized by a large enough field B + , along 
a desirable spin-quantization axis. This could pave the 
way for a gate- or strain-tunable spin filter in two dimen- 
sions, without the use of ferromagnetic materials, which 
may become useful for developing semiconductor-based 
spintronic applications. 



V. DISCUSSION 

Let us finally discuss some potential practical conse- 
quences of controlling the curvature field (24) by tun- 
ing the appropriate SO coupling profile A(rh The lo- 
cal SO strength A determines the curvature component 
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